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ABSTRACT 

In this paper, the RW Aur A microjet is studied from the point of view of X- 
wind models. The archived HST/STIS spectra of optical forbidden lines [0 1], [S 11], 
and [N 11] from RW Aur A, taken in Cycle 8 with seven parallel slits along the jet 
axis, spaced at (/.'07 apart, were analyzed. Images, position-velocity diagrams, and line 
ratios among the species were constructed, and compared with synthetic observations 
(■h I generated by selected solutions of the X-wind. Prominent features arising in a steady- 

Qh| state X-wind could be identified within the convolved images, full- widths at half maxima 

and high-velocity peaks on both of the redshifted and blueshifted jets. The well-known 
asymmetric velocity profiles of the opposite jets are built into the selected models. 
We discuss model selections within the existing uncertainties of stellar parameters and 
inclination angle of the system. In this framework, the mass-loss rates that were inferred 
to be decreasing along the jet axis in the literature are the results of slowly decreasing 
excitation conditions and electron density profiles. Despite the apparent asymmetry in 
terminal velocities, line intensities and mass-loss rates, the average linear momenta from 
the opposite sides of the jet are actually balanced. These previously hard-to-explain 
features of the asymmetric RW Aur A jet system now find a different but self-consistent 
interpretation within the X-wind framework. 

Subject headings: ISM: jets and outflows - ISM: Herbig-Haro objects -ISM: individual 
(RW Aur A) - stars: formation 



1. Introduction 

The RW Aur system is one of the brightest and mo st studied members of the Taurus- Auriga star 
forming region (V « 10.51: iKenvon fc Hartmaimlll995l ). It is a wide binary system with a separation 
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of ~ l'.'40 (jWhite &: Ghea l200lh . Both mem bers possess circumstellar disks which are tidally 
interacting with each other (jCabrit et al,l 120061 ) . The primary, RW Aur A, is the o nly source that 



drives a bipolar jet. T he jet was first inferred from its optical forbidden line profi les (IHamann 



1994 



Hartigan et al.l 119950 and later confirmed by narrow-band imaging (HH 229; iMundt &: Eisloffel 
19981 ). RW Aur A is known to be activel y accreting with an accretion rate of 0.5 — 1.6 x 10 -6 



M yr 1 inferred from optical veiling (e.g.. lHartigan et al. 



1995 



White fc Hillenbrand! 12004 ). 



Mundt fc Eisloffel 



Mundt fc Eisloffel 



Overall, RW Aur A jet has a strikingly well-collimated appearance (e.g., 

1998J: Dougados et al.l[200ol : lwoitas et alJboO^ . The [S II] AA6716+6731 image of 
(|199g ) showed that the blueshifted jet lies exactly along the position angle of 130° at least up to 110". 
The jet is marginally resolved within 56 AU by CFHT adaptive optics and HST/STIS. Widths of 



the jet slowly incre ase beyond this distance with an opening angle of 3° to 4° (IDougados et al 



2000 



Woitas et al.ll2002l ). and they appear to expand more rapidly closer to the source (within Cf.'5) than 
farth er away. Such trends of jet widths have been difficult to interpret in the past (jReipurth k, Bally 
200ll ) , as often the observed widths were interpreted as the extent of bow shock wings and the true 



jet should be much narrower. Farther than / /5 from the sources, the widths of the RW Aur A jet 
become wider as the in tensity decreases along the jet axis, which was also observed by HST for 
HH 34 (|Rav et al.lll996l ). The fluxes of both blues hifted and red s hifted jets decrease by an order 



of magnitude from Of! 5 to within 2" from the star. IWoitas et al 



widths margin ally consisten t with an X-wind jet in IShang et al 



wind model in lGarcia et al.l (|200ll ) at high accretion efficiencies. 



(2002) considered the measured 
(2002) and the self-similar disk 



Kinematics based on long-slit spectroscopy makes the RW Aur A jet a textbook example of a 
microjet. The position- velocity diagrams of its brightest forbidden lines, such as [S n], [O i], and [N 
n], clearly display one single intensity peak at each position along the jet axis. The centroid velocity 
lies around 100 kms -1 which is identified as the high-velocity gas tracing the collimated jet. The 
velocity widths are b road close to the source and narrow to its terminal v e locities far along the 
length of the jet (e.g 



Bacciotti et al 



1996 



Hirth et al 



1997 



Pvo et al. 



2006 



Melnikov et al 



2009 



Hartigan fc Hillenbrandlbood ). The low- velocity intensity peaks around 5 to 20 km s 1 (IHirth et al 



19971 ). which are usually presen t in many classical T-Tauri stars, are vir tually non-existent either 



in the optical HST/STIS data (IWoitas et al.ll2002l : iMelnikoy et alJl2009l; an d this work) or in the 
near-infrared [Fe n] A1.644 urn spectra fr om Subaru/IRCS (IPvo et alJboodl and Keck/N IRSPEC 
( Hartigan &: Hillenbrand! 2009 ) . Although Hartigan et al. ( 1995 ) and Alencar et al. ( 2005 ) reported 
detections of some low-velocity peaks in [O i] A6300, these peaks were highly variable and have 
been discarded as simply transient on a timescale of months. The strong collimated emission 
at high velocity makes the RW Aur A jet one of the ideal candidates to confront with models of 
magnetocentrifugally-driven winds arising from a n arrow region close t o the inner-edge of the disks. 



for ex ample, those that may resemble the X- winds (jShang et al.lll998l ; IShu et al.ll2000l ; IShang et al 



20071 ). In fact, it is our primary goal to compare this set of data with the best fits g enerated by 
models based on the X-wind solutions using approaches outlined in IShang et al.1 (120021 ). 



Despite the predominant high- velocity emission, asymmetry across the opposite sides of the RW 
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Aur A jet has long been documented in the l i terature (iHirth et al.lll994l. 119971: iBacciotti et al.lll996l; 

Woitas et allhood : iLopez-Martm et all20o3 : IPvo et al.ll2006l : lMelnikov et all2009l ; lHartigan k Hillenbrand 



20091 ) . The level of asymmetr y in RW Aur A is among the largest in the 8 asymmetric T-Tauri jets 
reported in lHirth et al.l ([19941 ). who excluded interpretations based on different spatial orientations 
on opposite sides of the source, uncorrelated time- variabilities, different degrees of deceleration, or 
different degrees of collimation. Instead, t hey concluded that such strong asymmetry should arise 
in the acceleration process of the jet itself. IWoitas et al.l (120021 ) derived a velocity ratio of 1.6— 1.76 
within Of! 5 from the source with the HST/ STIS, and attributed the origin of asymmetry to the 
"central engine" . ILopez-Martm et al.l (J2003J) measured proper motions of the knots and found that 
they have the same ratios as their radial velocities, and the velocity asymmetry should be intrinsic 
to the jet. The ratios of proper motion are maintaine d well up to 4f 6 and 1 .6 from the source, 
along the blueshifted and redshifted jets, respectively (jMcGroartv et al.l 120071 ). We will treat the 
velocity asymmetry as an intrinsic property built in our models to explain the overall phenomena. 

The asymmetry exists not only in flow velocities but also i n line intensities and physical 
conditions derived from the line ratios. iMundt fc Eisloffell (|1998l ) noted that the two opposite 
sides of the jet differ significantly in their spatial extent, with the [S n] jet extending to ~ 106" 
and ~ 50" o n the blueshifted and redshifted sides, respectively. The red jet is brighter in the 
inner ~ 10" (IBacciotti et al.lll996l ). and remain s so down to Of! 2 from the so u rce as seen in both 
optical (jMelnikov et alJbood ) and near-infrared (jHartigan & Hillenbrandlbood ) . Iffirth et al.l (|1997l ) 
noted that, at ~ 1" angular resolution, the red jet has a larger electron density, a higher [S n] 
A6731/[0 i] A6300 ratio, and a lower Ha/[S n] AA6716 + 6731 ra t io, th an the blue jet. On the other 
hand, the faintness of the blue jet prevented IBacciotti et al . ( 1996) from attempting a reliable 
diagnostic analysis using ground-based spectra, which led lWoitas et al.l (|2002l ) to assume similar [N 
II] A6583/[0 i] A6300 and [S II] A 6731/[Q il A630 ratio s for both sides of the jet for the HST/STIS 
spectra in their later analysis. Melnikov et al. (j2009l ) re-examined the same STIS dataset and 
obtained a smaller [N II] AA6548 + 6583/[0 i] AA6300 + 6363 but a larger [S II] AA6716 + 6731/[0 
i] AA6300 + 6363 in the red jet, which imply a lower ionization fraction and temperature than the 
blue jet. We will re-examine the dataset, and analyze the lin e ratios and physical conditions based 
on the diagnostic approach developed in IShang et al.l (|2002l ) for the X-wind jets. 



As a result, the fundamental physical properties of the RW Aur A system may appear asym- 
metric. Mass-loss and momentum rates will be asymmetric when they are inferred from the mea- 
sured asymmetric radial velocities, jet w i dths, and the electro n densities obta i ned f rom line ratios 
(|Bacciotti et al.lll996l : Iwoitas et alihood : iMelnikov et alibooa ). IWoitas et all faooi ) reported that 
if similar physical conditions had been adopted, the red jet would carry one half of the mass loss 
and one third of the momentum rates c arried in the b lue je t. The imbalance of momentum may 
cause dynamical effects, such as recoil. I Wait as et al.l ( 20021 ) calculated an upper bound 8kms _1 
for the recoil, under the assumption that the imbalance persisted during all the pre-main sequence 
lifetime of the star. They stated that this recoil value is probably a gross overestimate, because 
it is possible that earlier accretion of mass into the star might have been more symmetric, or that 
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linear momentum balance might have been maintained over time. The dynamical time of these jets 
(~ 10 5 years or less) is relatively short, diminishing the reasons to expect an imbalance p ersistent 
during the total of the pre-main sequence phase (about 10 6 years). iMelnikov et al.l ( 2009 ) revised 
the calculation by independently deriving electron densities and ionization fractions from each side, 
and derived similar mass- loss rates for the two sides. They attributed the asymmetry to inhom o- 
geneous ambient mat erial or to the magnetic configurations. Either work, IWoitas et al.1 (J2002J) or 
Melnikov et al.1 (|2009i ). leaves the momentum unbalanced in the presence of the observed velocity 
asymmetry. In this paper, we will undertake similar analyses based on the same HST/STIS dataset 
used in previous works, and demonstrate that the basic jet properties should be balanced in spite 
of the observed asymmetries in velocity and mass loss. 

The X-wind framework provides an attractive possibility for models confronting with high- 
angular resolution observations. Through the star-disk interaction, the stellar magnetic fields can 
truncate the disk at the co-rotation radius in steady state, where the field is trapped and pinched 
toward a small area called the X-region at the inner edge of the disk. The material from the disk 
can accrete onto the star via a funnel flow along the field lines connecting the inner disk edge 
and the star. Meanwhile, the excess angular momentum from the star can be transported to the 
X-region and material can be magnetocentrifugally ejected along the open field lines as the X-wind. 
Streamli nes of the wind flow fan out from the X-region, and gradually collimate along the axis of 
rotation. IShang et al.l (j 19981 ) demonstrated the kinematic features of an X-wind through synthetic 



spectra with broad line widths close to the base of t he jet which narrow to a profile of single velocity 
component upon approaching flow terminal speed. IShang et al.l (120021 ) self-consistently calculated 
the ionization and thermal conditions for semi-analytic X-wind solutions, and demonstrated how the 
excitation conditions can influence the synthetic images as "optical illusions" , finding a good match 
of line ratios from availabl e Herbig-Haro objects and microjets. Based on the works established 
in (IShang et al.lll998l . l2002l ). here we attempt a comparison of the X-wind models with those well- 
known features of the RW Aur A jet. 

In this paper, we demonstrate a step by step comparison of the high-angular resolution data 
with theoretical models generated within the X-wind solutions through the analysis of key observed 
features. RW Aur A is an excellent candidate for such purposes due to its characteristic single high- 
velocity feature that simplifies the modeling. In Section 2 the analysis of ffi>T/STIS archival data 
is addressed, and our parallel procedures for producing synthetic observable quantities from the 
models are presented. The observed properties of the jet extracted from the data: i.e., images, 
position-velocity diagrams, and the line ratios are discussed in §3.1. An illustrative fit to the 
RW Aur A jet is given in §3.2, and directly compared with the observational data. Finally the 
interpretation of these unusual features of RW Aur A are discussed in the theoretical framework 
based on the X-winds in Section 4. 
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Methods 



2.1. Archival HST /STIS Data 

We obtained HST/STIS spectra of RW Aur A from the Multimission Archive at Space Tele- 
sc ope (MAST). The spectra were taken on 2000 Decem ber 10 and have been previously reported 



ioolfte 



used with the G750M grating. The slit was displaced at a step of Of! 07 parallel to the jet axis at 
seven different positions, creating a total field of view of about Cf.'5 across the jet. Each exposure 
produced a two-dimensional spectral image of 52" long in the spatial axis and 6295 to 6867 A in 
the dispersion axis, covering the seven brightest emission lines: [O i] AA6300, 6363, [N n] AA6548, 
6583, [S n] AA6716, 6731, and Ha. The slit sampled the emission at a nominal rate of / .'05 pixel -1 
with an angular resolution of ~ O'.'l (FWHM), and a spectral sampling rate of 0.554 A pixel -1 with 



a 2.5-pixel line width, corresponding to 

A. 



25 km s per pixel and a FWHM ~ 65 kms -1 at ~ 6500 



The pipeline-processed spectra taken directly from MAST work reasonably well for our pur- 
poses. Hot pixels were present at specific positions on the CCD, but they fall primarily on line- free 
regions. Further reductions were performed to remove the baseline undulations caused by the 
wavelength rectification of undersampled stellar continuum and the contribution from the reflec- 
tion nebula. Each of the seven spectral images was first divided into three sub-images containing 
lines of [O i], Ha+[N n], and [S n], respectively. On each sub-image, we examined the rows contain- 
ing the redshifted jet extending up to 4" and blueshifted jet up to 2" from the star. For each row, 
we fit the baseline by a Legendre polynomial, up to the 10 th order, with IRAF/STSDAS3 GFIT1D 
and downhill simplex x 2 minimization. 

Position- velocity (PV) diagra ms were obtained f or each line and each slit, relative to the 
systemic velocity of +16.0 kms -1 ( Woitas et al.ll2002l ). and their transversely averaged PV dia- 
grams. [O i] A6300 emission from velocity higher than —235 kms -1 from the blueshifted jet was 
cut off due to wavelength coverage of the G750M grating. Line profiles were also extracted from 
the PV diagrams every / . / 05. Each profile was fit with a single Gaussian and a linear baseline in 
I RAF/STSDAS NGAUS SFIT with recursive x 2 minimization. On the intensity maps as constructed 
in lwoitas et al.l ( 2002 ). a single Gaussian was fit to the transverse profiles at every / .'05 with the 
Levenberg-Marquardt method for the amplitudes, intensity centroids, and full width at half max- 
imum. For each pixel of the 2D STIS stack, ratios between pairs of lines were used to construct 
line-ratio diagrams for [S II] AA6716/6731, [N II] A6583/[0 i] A6300, and [S n] A6731/[0 i] A6300. 
The data points were further arranged according to their spatial distribution. 



1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of 
Universities for Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation; 
STSDAS is a product of the Space Telescope Science Institute, which is operated by AURA for NASA. 
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In the following analysis performed on the model data cubes, operations and conditions applied 
in the observational analysis were adapted so that the model dataset can be properly compared 
with the observed dataset. 



2.2. Synthetic Maps 



Synthetic images and long-slit spectra are calculated for various X-wind configurations with 
methods and examples illustrated in Shang et al. (1998, henceforth SSG) and Shang et al. (2002, 
henceforth SGSL). We follow the approach adopted in SGSL to obtain a self-consistent thermal 
structure for each model, in which various ionization and recombination and heating and cooling 
processes associated with the wind are considered. These include hydrogen photoionization and H 
photodetachment by photons from the star and accretion hot spots, X-ray heating and ionization, 
collisional ionization, ambipolar diffusion heating, and external heating arising from mechanical 
origins in the flows. In a magnetized star-disk system, the closed stellar magnetic loops confine 
plasma at soft X-ray temperature, and sites of magnetic field line reversals between the star and 
the wind are prone to reconnection events that can generate harder X-rays (see, e.g., Figure 1 of 
Shu et al. 1997). SGSL concluded that X-rays are effective in ionizing the base of the X-wind 
(see also Shang et al. 2010, henceforth SGLL), and mechanical heating can easily heat the flow to 
8000-10000 K with a small amount of mechanical energy converted to heat. In this work, we adopt 
conventions defined in SSG and SGSL, and find solutions normalized by stellar parameters for RW 
Aur A system. 

We extend the reaction network to include the ionization balance of oxygen, sulfur, and nitro- 
gen. Charge exchange with the H atom controls the equilibria between O and + , and N° and N + 
as their ionization potentials are similar (13.60, 13.61, and 14.53 eV for H, O, and N, respectively). 
Sulfu r is assum e d to b e ionized as a result of its lower ionization potential (10.36 eV). Atomic data 
from 



Mcndoza ( 



19831) were used in these calculations, and updated collisional strengths of [S II] 



(|Cai Pradhanlll993l ) were adopted 



To ob tain optimal fittin g with the observed spectra, the parameters for mass-loss rate M w 
defined in IShu et al.1 (|1994al ). X-ray ionization Lx (Equation 3.10 in SGSL), and coefficient of 
heating ah (Equation 5.2 in SGSL) were adjusted to obtain proper thermal profiles for emission. 
The fiducial cases shown in SGSL and Shang et al. (2004, henceforth SLGS) were taken as initial 
model standards. The initial values of M w and Lx were first allowed to vary with the ratio 
L x /Af w = 3 x 10 13 ergg" 1 fixed, a measure of X-ray attenuation previously in SGSL and SGLL. 
The values of are adjusted to produce thermal profiles that would generate emissions fluxes 
matching the observed ones. The properties of the synthetic images are derived and spectra for 
further analysis and comparisons with the STIS data. We note that the actual stellar system always 
looks knottier in emission than the synthetic profiles produced by the steady-state models, and that 



the level of X-rays adopted in our m odels in this pa 



mass pre-main sequence stars (e.g.. iGetman et al 



jer is on the higher end of flare s attained by low- 

201 ll ). As explained in 



2008 



McClearv fc Wolk 
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SGLL, the higher level of X-ray luminosity adopted in our modeling scheme may reflect the fact that 
these X-rays may have experienced higher extinction near the base of the wind to reach the level 
of ionization obtainable for an otherwise lighter wind. A lower X-ray luminosity can be adopted 
instead if the X-rays are generated at a more elevated position relative to where the ionization 
actually takes place. This can effectively lower the absorption column in the context of X-wind 
models (unpublished work). By adopting and keeping the fixed ratio Lx/M w = 3 x 10 13 ergg -1 , 
the analysis shown in this work is consistent and can be easily compared with past works using this 
approach (SGSL, SLGS, and SGLL). We define and focus on the overall trends and profiles rather 
than local variations within individual knots in the process of model-data comparison and fitting 
process. 

The synthetic images have to be generated at very high-resolution using high sampling rates 
for convolution. The maps were convolved with the telescope point spread function while retaining 
spatial sampling in the original map. STIS instrumental profiles were applied to the long-slit 
spectra. Each of the profiles was convolved with a two-dimensional Gaussian of Cf.'l x 65 kms -1 
FWHM, and resampled at Of! 05 and 25 kms -1 , based on profiles of STIS 52" x Of.'l slit and G750M 
grating. The resulting PV diagrams were fitted with a one-dimensional Gaussian at each spatial 
sampled point using the downhill simplex method. The fitted values of line intensities, velocity 
centroids, and velocity widths were compared with those derived from the observed spectra along 
the jet axis. 

The synthetic line ratios of [S n] AA6716/6731, [S n] A6731/[0 i] A6300, and [N II] A6583/[0 
i] A6300 were calculated diagnosing the jet from different positions in the integrated line intensity 
maps. 



2.3. Properties of the RW Aur A System 



Properties of the RW Aur A system affect the actual model fits to observations as the stellar 
parameters determine the actual dimensional values of the models. The systematic velocity also 
affects the values of the observed jet speeds relative to the driving source and the velocity ratios. 
Furthermore, the inclination angle may introduce uncertainties to the terminal velocities. Stellar 
parameters in the literature determined from various evolutionary tracks, veiling r\ and visual 
extinctio n Ay, stellar mass M*, stellar radii R*, spectral types and luminosity are summarized i n 



extinctio n Ay, stellar mass M t , stellar radii it*, spectral types and luminosity are summarized n 
Table fTI dHartigan et al.lll995l : Istout-Batalha et alJlinonl : Iwoitas et al1l200ll : Iwhite fc Ghe^boOlh 



It can be seen that the ranges of inferred stellar masses, radii and luminosities have not varied 
more than a factor of two. Within the uncertainties, we adopt a reference value of = 1M@ and 
-R* = 2i?Q for illustrative purposes throughout this work. 

RW Aur A m ay have a varying heliocentric radi al velocity. For th e same STIS dataset, 
Woitas et al.l (|2002l ) adopted a value of +16.0 kms -1 (IPetrov et al.l l200lh and obtained an av 



erage velocity ratio for the jet components of ~ 1.6. On the other hand, iMelnikov et al.l (120091 ) 
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adopted +23.5 k ms" 1 (Woitas et ali 2005 ) and obtained a higher value of the ratio ~ 1.8. The 
value obtained bv lWoitas et al.l (|2005h was based on one single fitting to the Li I A6708 absorption 
in the current HST/STIS dataset, which we reanalyze with a value of +22.1 + 0. 61 kms" 1 . For the 
illustrative purpose of this work, we adopt the averag ed value of +16.0 kms" 1 (IPetrov et al.ll200ll 



White &: Hillenbrandll2004l ) and set the nominal velocity ratio to be 1.6 throughout this work. The 
radial velocities reported in the literature are summarized and their inferred jet velocity ratios in 
Table [2j 

Various group s have attempted to derive the inclination angle of the RW Aur A system. 



(1998) and HST/STIS (2000) observations and derived a value of 53°. lLopez-Martm et al 



Looez-Martm et al. 


size. 


Alencar et al. 



(2003) 



2005|) 



modeled the emission and absorpt ion profiles of Balmer lines and Na I D lin es with a dis k wind 



model of lBlandford &l Pavnd (|1982l ) and inferred a range between 55° and 65°. ICabrit et al .1 ((2006) 



modeled 12 CO (J = 2 — 1) spectra of the disk around RW Aur A with Keplerian disk profiles 



following iBeckwith fe Sargentl (jl993l ). and found that the line profiles can be matched with 45° or 
60°. The range of inclination angles agrees reasonably well within 45° and 60°, which we will take 
as the uncertainty introduced when one de-projects velocities. 

Within the range of stellar parameters , we explore com binations of wind mass- loss rate M, 
disk truncation ra di us (the X-point, s ee 



( Shu et al 



1994 



30 



1995 



Shu &: Shane 



199 



Shu et al 

IT 



1994aT ) , and solutions of the X-wind models 



Shand 119981 : iShang et al. 



1998 



2002, 



2004 



2010 



[SSG, SGSL, SLGS, SGLL]: ICai et al.ll2008l ). Details of the magnetic configurations and their 
observational properties will be presented in a separate publication. In this work, we focus on 
specific features that one can derive from the observed properties of the RW Aur A jet, which are 
applied as selection criteria for the models. The X-winds that satisfy the identified constraints are 
further analyzed and compared with the STIS dataset. 

The velocity asymmetry, if intrinsic to the system, reveals an interesting connection of winds 
launched from the opposite sides of the disk in the X-wind context. In steady state, the inner 
disk truncation radius coincides with the X-point, and is co-rotating with the star at = J7 X = 

1 /2 

(GMx/R^) . T his intro duces the unit of velocity, v x = R x fl x , defined at the base of the jet 
dShu et al.l Il994al lbl Il995h . where the velocity of the wind can be non-dimensionalized as v w = 
Yb,rv x - The velocity ratio, v Wt B/v w ,R ~ 1-6, suggests the launch of unequal winds above and 
below the disk from the X-point. The R x is set as one single uniform value for both sides of the 
disk around the circumference. For the scope of this paper, we do not consider detailed launch 
mechanisms that could introduce the asymmetric mass loading or dynamical effects which may 
warp the disk or cause a non-circular orbit at the innermost edge. Such an investigation shall be 
reserved for future studies. 



We select X-wind solutions whose properties satisfy these conditions from IShangl (|1998l ) . The 
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ratio of the averaged velocity between the blue and the red jet components yields 

Vw,R T fi v x t R 

Without specifying the actual value of i? x , the selection of solutions will not be unique. In Tables 
and HI we give two sets of examples. The range of R x can vary from a value smaller than 0.05 
AU, to as large as 0.2 AU. Note that we have set no a priori condition that restricts the choices of 
R x by either stellar or model parameters. The reality probably exists somewhere between the two 
extreme values. Although further constraints may be available from other clues of the system, we 
leave it subject in our followup investigations. 

In Section 3, we analyze and compare with the archived HST/STIS data based on one choice 
of R x , one inclination angle, and one pair of solutions. The blue and the red sides are treated 
independently, although the two sides of the jet are connected through one common point R x in 
space and the fixed ratio of averaged velocities. The equiv alent indiv i dual m ass-loss rates M w 
defined in the sense of the original X-wind parameters (e.g., Shu et al. 1994al lbl) are adjusted to 



match the average intensities of knots in the respective side. The actual mass-loss rates that enter 
into either hemisphere are half of the inferred model values. The individual mass-loss rates for each 
hemisphere are reported in § 3 and § 4. 



3. Results 

In this section, an example to demonstrate the general principles in applying the X-wind models 
to observations is illustrated. We first show the observed features of the RW Aur A jet extracted 
from the archival HST data. 



3.1. Observed Features 

The two sides of the RW Aur A jet are asymmetric in appearance. The narrow-band images 
of three forbidden lines [S n] A6731, [N II] A6583, and [O i] (A6300 for the red jet and A6363 for 
the blue jet) are shown with gray scale in Figured! Both the blue and red jets are knotty, and the 
knots become fainter as the jet propagates outwards. The innermost knots appear to be tightly 
connected while the distance between the outer knots increases. The red jet extends beyond 4" 
and at least four knots can be identified. The blue jet is overall dimmer by factors of 3 to 4, and 
extends up to 2" with three knots identified. The extent of asymmetry slightly differs among line 
species. [N n] is overall dimmer and shows less asymmetry than [O i] and [S il], and lacks emission 
beyond R4. The red jet appears somewhat more compact in the transverse direction than the blue 
jet. 

The asymmetry in intensities and velocity centroids are indeed pronounced features on the PV 
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diagrams. The velocity ratios between the blue and red jets are ~ 1.56 to 1.60, which can be traced 
to the innermost detectable region of the system at 072. Figure [2] shows the PV diagrams of [O i] 
(A6300 for the red jet and A6363 for the blue jet), [N II] A6583, and [S II] A6731 extracted from 
the central slit. The velocity ranges from —25 to +300 kms -1 in the red, and from —350 to +25 
kms -1 in the blue. The knots were identified by peak positions in the PV diagrams and are labeled 
as in Figure [TJ The line shapes of the knots can be fit with a single Gaussian profile throughout 
the jet, showing clearly that the jet is dominated by a single velocity component. The velocities at 
the centroids range between —160 and —180 kms -1 and +100 and +120 kms -1 , on the blueshifted 
and redshifted sides, respectively. Velocity centroids within the observed field of view vary slightly 
among the species. The velocity along the jet varies, as the knot R5 in the red jet is ~ 10 kms" 1 
faster than the overall average, and the knot B3 in the blue jet shows a decreasing trend on the 
order of ~ 20 kms -1 . 

Loci of line ratios from the blue and red jets do not overlap on the cross-ratio plots. The 
line ratios between [N n] A6583/[0 i] A6300 and [S n] A6731/[0 i] A6300, and between [S n] 
AA6716/6731 and [S n] A6731/[0 i] A6300 are shown in Figure El We show the line ratio plots for 
points within / . / 95 (B4+B5) from the source for the blue jet, and within / .'65 (R7) for the red jet, 
at a distance ratio approximately proportional to their respective velocities. Line ratios integrated 
over slits labelled by their positions from the source are shown in the left panels, and those from the 
individual slits are shown in the right panels. Only those points with relative uncertainties lower 
than 50% are shown. We focus on the analysis for data points from the three central slits since 
the outer slits do not contribute statistically significant information as seen in the right panel of 
Figure [31 The line ratios increase roughly with distance along the jet axis. The line ratios derived 
from different slits do not segregate into distinct groups, but rather scatter around a general trend. 
The emissions appear to come from continuous media of similar physical conditions within each 
side of the jet, but each component has its own distinct trends. In the red jet component, the 
[S n] AA6716/6731 ratio clusters around similar values 0.45 — 0.63 whereas the [S n] A6731/[0 i] 
A6300 ratio ranges from 0.3 to 1.0. It also shows an increasing trend with increasing values of [S 
n]A6731/[0 i]A6300 from the outer slits. In the blue jet component, both the ratios of [S n] doublets 
and [S n] A6731 to [O i] A6300 have large scatters due to the weakness of the lines, and there is no 
clear correlation of the values. Similar trends exist for the [N n] A6583/[0 i] A6300 ratio: those in 
the blue span a wider range while those in the red are located within a confined region, however, 
the two line ratios are clearly positively correlated, from the central to outer slits. Perhaps due to 
the faster drop-off of [O i] signal particularly in the blue jet, points from farther distance tend to 
drift and pile upward in the upper-right corner of the [N n] A6583/[0 i] A6300 to [S n] A6731/[0 
i] A6300 panel with even larger scatter. Those points were not displayed in Figure [3] for clarity of 
presentation. 
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3.2. A Sample Fit for the RW Aur A System 

In this subsection, a sample fit for the RW Aur A jet is shown. For simplicity, we demonstrate 
one case that could capture overall the major features of the jet in this paper. Other solutions are 
possible, within the observational constraints set by the HS T/STIS data, such as those listed in 
Table H 

The sample is described by properties outlined in Section 2.3 for the RW Aur A: stellar param- 
eters (M*,/?*) = (1M0,2J?q), a systemic velocity of +16.0 kms" 1 (Table [2]), an inclination angle 
of 55° (Table [3]), and a disk truncation radius i? x = 0.12 AU (see Table [3] for the range of possible 
choices). With these parameters, one solution gives an average projected velocity of —170 kms" 1 
for the blue side and the other gives +115 kms -1 for the red side as inferred from the archival 
data. The velocity coverage is in the range of 200 to 310 kms -1 for the blue, and 120 and 200 
kms -1 for the red. They possess properties of those pointed out by SSG such as the cylindrically 
stratified density structures, broad line profiles near the bases, and slowly-collimating streamlines 
toward the axis over a large distance. The density contours and streamlines are similar to Figure 
1, and the PV diagrams are similar to Figure 3 of SSG. 

We follow SGSL in finding the optimal parameters for mass-loss rates M w , X-ray luminosities 
Lx and mechanical heating coefficient a^. During the parameter searches, the ratio of X-ray lumi- 
nosity to mass-loss rate, Lx/M w is set to be 3 x 10 13 erg g -1 , as in SGSL, SLGS, and SGLL. From 
parameters adopted in the fiducial case of SGSL for both winds in the blue and red hemispheres, 
M w , Lx, and at are individually varied on each side to find the optimal matches for the intensity 
profiles and line ratios. The optimal sets of parameters are given in Table which capture the 
overall features reasonably well. Although the emission is knotty and the profiles and velocities 
have variations and scatter, clear average trends can be clearly identified for individual lines along 
the jet axis. For individual mass-loss rates of 3 x 10 _8 Mq yr" 1 in blue and 5 x 10" 8 Mq yr -1 in red, 
and ah = 0.0010 for the blue and = 0.0015 for the red, the physical properties and emission 
profiles are illustrated in Figures H] to 

The physical scale in Figures [H and [5] covers the transverse and longitudinal dimensions of the 
seven slits: 50 AU in horizontal distance w and 150 AU in length z along the jet axis. Figured] shows 
the velocity contours on top of the electron density profiles, and Figure [5] shows the distributions of 
electron fraction and temperature. In Figured! the electron density is also cylindrically stratified 
as the underlying density as described in SSG from within 50 AU out to larger distances. Near 
the base of the jet, it casts an opening of a conic shape in the images as pointed out in SGSL. 
The velocity increases toward the axes on b oth sides of the jet, where the flows are denser and 



faster, until reaching the terminal velocities (jShu et al.lll995l : SSG). In Figure [21 electron fraction 



and temperature follow similar behaviors to those in SGSL. Both x e and T increase toward the 
inner boundaries of the wind, and peak within 5 AU from the launch point. The blue jet has a 
higher level of electron fraction and temperature than the red jet, supporting the trends observed 
from the line ratio diagnosis diagrams. 



-12- 



We show comparisons in PV diagrams, intensity profiles, image full-widths at half maximum 
(FWHM), and velocity centroids in [0 ijj [N II] A6583, and [S n] A6731 in Figures © to [SJ On 
the left, the profiles of intensity peaks and their FWHMs along the jet are shown for the inner 2'.'5 
from the source. On the right, the contours of the synthetic emission spectra (red for the red jet, 
and blue for the blue jet), convolved with the instrumentation profile, are overlaid on top of the 
observed PV diagrams. The lowest contours and contour intervals of both spectra are identical to 
those of Figure El Velocity centroids from both the observed and synthetic spectra are shown at 
the STIS velocity resolution of ~ 65 kms -1 . The averaged velocities at centroids of the synthetic 
spectra are ~ —170 kms -1 for the blue and ~ +110 kms -1 for the red, which reconfirm the velocity 
ratio of ~ 1.6, and the trend of slightly increasing velocities along the axes. In the case of [N II] 
from the blue jet, emission at the stellar position is missing. However, we can not assess whether 
this is due to the nature of weaker emission from the blue jet or other reasons. 

The line intensities and image widths of the synthetic images cover the overall trends in the 
observed ones. The trends delineated by the theoretical curves capture the major features well up 
to locations where the observed knots become too faint. For the redshifted jet, the model profiles 
go through the average of the knot peaks and off-knot regions; in the blueshifted jet, the same 
holds for the intensity profiles, although for the image widths and velocity centroids, the curves 
go through more scattered points except for the [S Ii] A6731. The image widths increase gradually 
with distance as the overall intensity profiles decrease gradually along the axes. This is expected 



as the electron fraction overall follows a pattern of recombination pointed out by iBacciotti et al 



4l996j) and SGSL. 



The similarities and differences in the cross line-ratio diagrams between the synthetic maps 
and the actual data reveal the degree to which the physical conditions are reproduced by the model 
parameters. Figure [9] shows the line ratios between [N Ii] A6583/[0 i] A6300 and [S n] A6731/[0 
i] A6300, and between [S Ii] AA6716/6731 and [S Ii] A6731/[0 i] A6300 overlaid by model coverage 
(in color shades). The shaded areas include the innermost Cf.'95 of the blue jet and 0'/65 of the red 
jet, at a distance ratio approximately proportional to their respective velocities. In the transverse 
direction, the model points extend to about w = / .'l2, corresponding roughly to the spatial coverage 
of the central three slits. The shaded areas generated by the synthetic images encompass most of 
the observed data points (in discrete symbols), following similar trends, up to the tips of the knots, 
B4 and R7. The overall ratio of [N n] A6583/[0 i] A6300 is about twice as large in the blue jet as 
in the red jet. As shown in Figure [5j the overall electron fraction of the blue wind is approximately 
twice as large as in the red wind. The spatial dependence of [S n] A6731/[0 i] A6300 varies with the 
differences in critical densities in [S n] A6731 and in [O i] A6300, and the [S n] doublet ratio varies 
across the jet as n e decreases faster laterally than longitudinally in the model. The combination of 
(x e ,T) in Figure [5] following our selection criteria of the solutions indeed captures the overall sense 



2 We remind the reader that different [O i] doublet lines are used for the red and blue jets: A6300 for the red and 
A6363 for the blue, due to the cut off of the blue wing of A6300 at instrument level. 
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of differences and asymmetry found on the opposite sides of the jet. The differences in excitation 
conditions within the opposite jet components differentiate the distribution of data points on the 
respective line ratio diagrams. 

Our synthetic points are distributed over a wider area than the observed points. This is 
because a large area is covered by fainter pixels on the computational domain. This is especially 
noticeable for the brighter red jet. The model points come from the equivalent area covered by 
the central three slits as small pixels. Emission from these theoretical pixels is excited ideally by 
the model approach outlined in § 2.2. As each of the theoretical pixels contributes to the model 
line-ratio points, there may not be corresponding points from the observed data. This is because 
the actual jets are knottier than the model ones and some of the jet volume in real space is not 
sufficiently excited as in the model counterpart. This explains why the model naturally predicts 
a broader range of values in line ratios while the observed range is narrower because observations 
are naturally weighted toward the regions with strong emission. On the other hand, the observed 
data for the blue jet have much higher uncertainties because the jet itself is weaker in emission 
but faster in flow speed, which may result in a wider spread in physical conditions to begin with. 
However, a rough trend can still be traced as in the left panel of Figure El when all the slits are 
added to increase the signal-to-noise ratio. In this case, an overall match in the global trends are 
what we are seeking in the process. We also note here that model points from the equivalent of 
three central slits have already covered most of the regions where emission lines are excited. Adding 
the emitting area of two more slits outside of the central area does not produce significant increase 
of the emission. Rather, the additional emitting area will mostly contribute to points with lower [O 
i] A6300 emission, resulting in higher [S II] A6731/[0 i] A6300 and [N II] A6583/[0 i] A6300 ratios, 
and [S n] doublet ratios approaching unity. These points may not be sensitively sampled compared 
to the brighter points represented by the observed data. 



4. Discussion 



4.1. General Features 



We hav e attempted to re-a nalyz e the set of archived ff £T/STIS data previously analyzed and 
published in I Wait as et al.l (|2002l ) and lMelnikov et al.l (|2009l ) . Similar results were reached in general, 
although with some noticeable differences in terms of infer r ed ph ysical quantities. Our methods of 
data analys is differ from t hose ad opted in Melni kov et al.1 (|2009l ) , but are mostly similar to those 



Melnikov et al 



!00j|) included further noise subtraction and 



adopted in IWoitas et al.1 (|2002l ) . 
correction of Ay ~ 0.4 throughout the jet. The remaining noise in our data could introduce more 
scatter for the line ratios in the blue jet, but the overall values are co nsistent since the line ratio s 
in adjacent wavelengths can mitigate the effect of extinction to ~ 10% (jBacciotti fc Eisloffellll999l ). 
They also adopted a higher radial velocity for RW Aur A (+23.5 kms -1 ) than the averaged value 
(+16.0 kms -1 ). As expected, the line properties are nearly identical apart from the systematic 
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velocity shift. On the results from archived data, the overall properties one would infer in general 
from the opposite sides of the jet are not altered by different approaches and assumptions. 

The physical conditions derived from the different diagnostic analyses, r ather than the method s 
of data analysis, distinguish the various works based on the sam e set of data. iMelnikov et al.l ([2009) 



derived the physical quantities using techniques developed by iBacciotti fc Eisloffe 



1999) for the 



individual jet w i thout making a priori assumptions of the fainter blue jet as in lWoitas et al.l (|2002l ). 
Melnikov et al.l (|2009i ) reached the conclusion that the blue jet is less dense and more excited than 
the red jet. The abundances of the atomic and ionic species involved are somewhat different, 
which may contribute to theoretical line intensities and ratios through the combination of electron 
density, ionization fraction, and temperature. Comparing the physical conditions derived from 
our self-consistent approaches (Figures [4| and El) and those d erived from diagnostic mapping by 
employing the BE technique (Table 1 of IMelnikov et al.ll2009l ). we note differences in values and 
relative trends. For example, electron densities derived from the BE technique are lower in the red 
jet by 20%. The faintness of the blue jet introduces larger uncertainties in line ratios and in the 
derived electron densities. Values derived from the [S n] doublet may also represent a lower limit 
due to its low critical density (~ 2.5 x 10 3 cm" 3 ). It can be seen in FigureUJthat values for both jets 
are similar in the densest region, but they decrease faster in the cylindrical distance w in the blue. 
From the ionization fraction and temperature, both works concluded that the blue jet is the more 
heated and ionized, and the temperatures derived from the BE technique are typically 50% higher 
than those found in this work. The higher temperature may compensate for the lower inferred 
electron densities for the observed line intensities. The re- heating of the jet close to knot positions 
may raise the values locally. The derived ionization fraction from IMelnikov et al.l (120091 ) is in line 
with our values at the innermost streamlines where X-ray irradiation is the strongest (Figure [5]) . 
The hydrogen density that is derived by dividing the electron density with the electron fraction, 
may weight toward the lower value, leading to mass-loss rates that could be lower by an order of 
magnitude (see below). 

Physical properties such as the mass- loss rates inferred for the opposite sides of the jet are 
the best tests for the dynamical models and diagnostic analyses developed and in tegrated for the 
model s. Our best-fit mass- loss rate for t he red jet, is in line w ith t hat found in 



( 19961 ) using line ratios, and those found in lHartigan et al.l ( 19951 ) and lWhite fc Hillenbrand! ( 2004 ) 



Bacciotti 



using forbidden line intensities. These values are of the sa me order as thos e in 



Woitas et al 



et a 



(200 



2 



bu t one order of magnitu de larger than those derived by Melnikov et al. ( 20091 ). The small value 
Melnikov et al.l (|2009l ) lies in their consideration of "intrinsic" image widths. Although the 



jet radii appear more accurate, the mass loss carried by the associated wide-angle wind is not 
taken into account. It is also possible that the electron densities derived from [S n] doublet may 
represent a lower limit. On the other hand, the imbalance of momentum de d uced from the ratio 
of mass-loss r ates b etween the two sides of the jet is larger in IWoitas et al.l (J2002J). Specifically, 



Woitas et al.l (J2002J) obtained line-of-sight velocity, vj, and jet width, rj, from the same dataset 



and derived electron densities n e from the [S n] line ratio. They inferred a mass-loss rate M w 
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MfftH ( n e I Xe) r KVjT 1 ; I cos i, where electron fraction x e of the red jet was adopted from lDougados et al 



(|2002l ). and the values in the blue jet were assumed to be similar to those in the red jet. They 
also assumed similar [N n] A6583/[0 i] A6300 ratios for both sides, which led to an estimate of 
a mass- loss rate ratio ~ 2 (blue to red). Similarly, a ratio of ~ 3 is derived for the momentum. 



Woitas et al.l (|2002l ) interpreted this as a temporary phenomenon during accretion since there was 
no recoil velocity detected. Yet the discrepancy may be resolved if the adopted x e is 2 to 3 times 
higher in the blue jet, then the ratio of momentum rates from the two sides would be, oc v? /x e ~ 1, 
after correcting the ratio of electron fraction. In addition, their analysis of the emitted fluxes 
gave mass- loss rates that decrease by more than one order of magnitude within a distance of 2", 
which they interpreted as a decrease in mass loss due to sideway push at internal working surfaces. 
Although they also acknowledged the possibility of a decrease in electron fraction and a reduction 
of excitation, they left the unbalanced mass-loss and momentum rates unsolved. Our method self- 
consistently calculates the physical properties that best match the overall line intensities and line 
ratios, not biased toward only the inner dense part, leading to a better estimate of the mass- loss 
rates and a general balance of linear momentum. 

The "excess" redshifted emission at the stellar positi on on the approa chin g side is difficult to as- 



sess d ue to subtraction artifacts of the stellar continuum. |Pvo et al .1 (|2006l ) and lHartigan fe Hillenbrand 



(|2009l ) pointed out the high- velocity features of the RW Aur jet in [Fe n] A1.644 /jm at a velocity res- 
olution < 30 kms -1 : narrow line widths at the jet proper, symmetric profiles at each position, and 
increasing line widths toward the source, in general supports the predicted PV diagrams shown in 
SSG, rather than the disk-wind models. On the other hand, they also commented that the original 
PV diagrams for an approaching X-wind jet shown in SSG perhaps carry some "excess" re dshifted 



emission extending over the stellar position toward the counter side. lHartigan fe Hillenbrand! (120091 ) 
suggested stronger heating (or re-heating) at larger distance may resolve the discrepancy. The 
wide-angle wind would inevitably possess streamlines that contribute to the redshifted line-of-sight 
velocity, if the system is not strictly pole-on. We followed SSG and SGSL for the calculations of the 
line emissivities. The so-called "excess" is less pronounced in our maps generated for RW Aur A 
parameters with self-consistent excitation conditions implemented, when compared to the original 
PV diagrams in SSG that were produced with uniform temperature and ionization throughout the 
wind. This might suggest that the inner streamlines in our system might have been slightly more 
excited than the wind case inside the RW Aur A system. With O'.'l spatial resolution of HST, the 
emission is restricted to the stellar position within one resolution element, and this area is strongly 
affected by the artifa ct introduced by stel lar baseline subtraction in a rectified two-dimensional 
STIS spectra (see also iMelnikov et al.ll2009l ). It is very difficult to assess the actual "excess" of the 
redshifted emission. 

We have demonstrated the existence of a pair of X-wind solutions that could capture the 
major features of asymmetry in the RW Aur A system. The asymmetric properties reported in 
the literature and the archival HST/STIS data are confirmed by our analyses. Mass- loss rates 
that are required to best match the overall trends in the emission profiles and line ratios are also 
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asymmetric: Mb = 3 x 10~ 8 M yr _1 for the blue, and Mr = 5 x 10~ 8 M & yr _1 for the red. These 
two mass-loss rates give a ratio of Mb /Mr ~ 0.6. Multiplied by the asymmetric terminal velocity 
ratio, vb/vr ~ 1.6, we reach a ratio of linear momenta Mb vb/Mrvr ~ 1. This suggests that 
momentum ejected into the opposite hemispheres is indeed balanced. This is in contrast to the 
apparent d ilemma of unbalanc ed momenta carried by th e two jets from the velocity asymmetry 
reported in lwoitas et al. ( 2002 ) and Melnikov et al. ( 2009 ). Despite the asymmetry in velocity, line 
intensities, and mass-loss rates, the system is on average balanced for the most important physical 
quantity. 

The apparent unequal mass-loading and the uneven terminal velocities in the northern and 
southern hemispheres of RW Aur A do not break the law of conservation for linear momenta 
and magnetic fluxes contributed in the wind in our study. These major physical quantities are 
indeed conserved in the average sense across the opposite jets. Prom the analysis of emission lines 
originating from the jet alone, the conservation is an integral pa rt of the origina l s cenari o and 
framework laid out in the original works of the X-wind models ijshu et al.1 [l~994aU bl ll99fJ ). and 
there is no need to invoke asymmetry beyond the immediate wind proper. There is also no obvious 
detectable precession or signatures from precession of the jet that would suggest a non-alignment 
of the magnetic axis with the rotation axis on the physical and time scales of the existing datasets. 
However, we can not rule out other possibilities from the star itself or processes of the star-disk 
interaction. On the other hand, if the timescales of variation are short, effects of asymmetry 
originating deep in the star system might have been averaged over the rotation periods and long 
flow time so the visual effects on the large scale are small. Further clues along these lines supported 
by other observational evidence will be discussed in a follow-up investigation. 



4.2. Further Constraints 

In this work, we demonstrated the possibilities of finding a match between the X-wind models 
and the HST/STIS spectra of the RW Aur A jet, and discussed additional properties that were 
inferred from such a matching process. Although the existence of such a set of solutions does not 
constrain a unique launch point R x , the overall procedure suggests that general features of the X- 
wind models indeed match those obtained from the data analysis. In particular, the mass-loss rates 
and excitation properties not only fall nicely within the general range suggested by the data but also 
solve the puzzling asymmetry in the physical properties of the system. Despite such ambiguity, the 
general results demonstrated in this work should exist in any model that is to be selected should 
a different combination using another i? x provides a good fit with the detailed properties. The 
differences should be within a factor of two of what we learned in this work. 

Current available observational data on the RW Aur A jet alone do not seem to be able to 
resolve the ambiguity. Observational clues or constraints derived near the base of the jet related 
to the inner truncation radii of the disk, or some other information inferred from the permitted 
lines originating in the accretion funnels may help resolve the uncertainties. Inner truncation radii 
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of the RW Aur A disk from near-infrared interferometry (jAkeson et al.l 120051 ; Eisner et al.l 120071 ) 
have a range of 0.08 to 0.2 AU, consistent with fits of models with d ust sublimation tem perature 
between 1000 to 1500 K. The gas disk may extend to 0.02 to 0.05 AU (|Eisner et al.ll2009l ) inside of 
the dust disk although the actual values may be dependent on the underlying disk models adopted. 
These are all well within the considerations of choices for f? x given in this paper (Tables [3] and |4|) . 

Recent Zeeman-splitting measurement of optical emission lines can help to qua litatively in 



fer th e magnetic structure of the interaction region between the star and the disk (jDodin et al 



2O120 . The averaged longitudinal magnetic field strength B z inferred from the narrow component 
of He I A5876 was found to change pola rity within one rota tion period assuming P = 5.576 days. It 
is consistent with a model proposed in iPetrov et al.l (|200ll ) that the variation of the lines is due to 
two hot spots on opposite poles of the star. The co-rotation radius corresponding to this period is 
~ 0.07 AU assuming M* = Mq. While this is still within our considerations, it may also be possible 
that the disk truncation radius ( ~ R^) does not always coincide with the co-rotation radius. In the 



picture of a fluctuating X-wind (|Shu et al.lll997l ). f2* 7^ il x , the magnetic energy released would be 
proportional to the difference in angular frequency, and part of it can be in X-rays. In this case, 
we know that the X-ray activities may affect the initial ionization of the flow and properties of the 
lines. Simultaneous monitoring in both the X-rays and photospheric and forbidden emission will 
be desirable to understand the star-disk interaction in the RW Aur A system. 



5. Summary 

We have presented detailed analyses and comparisons between the HST/STIS archival data 
and synthetic spectra from the X-wind models for the bipolar jet of RW Aur A. We analyzed 
properties of optical forbidden lines [O 1], [N 11], and [S 11] through PV diagrams, images, and line 
ratios. PV diagrams of all the lines show similar kinematic features, and an average blue-to-red 
velocity ratio of ~ 1.6. Spectra of high spatial resolution enable diagnostics for both jets through 
the line ratios. We were able to identify general properties related to the asymmetry of physical 
conditions derived through the line ratios, and applied them as constraints for model comparisons. 

The diagnostic approach developed in SGSL was extended and applied to RW Aur A, using 
parameters for the stellar mass, radius, systemic velocity, and inclination angles of the jet and disk. 
Within observational uncertainties, we could select a set of stellar parameters to demonstrate a 
sample fitting process for illustration. Parameters were searched for mass-los s rates, X-ray lumi 



nosit i es, an d mechanical heating starting from the fiducial models adopted in IShang et al.l (|2002 



20041 . |2010| ) for a pair of X-wind solutions that were selected by matching the velocity ratio. The 
combined models not only possess the feature of asymmetry in jet speeds, but also simultaneously 
meet the overall agreement in line intensities, image widths, and line ratios. 

Our interpretation of the RW Aur A system offers a different perspective to the driving of 
the asymmetric jets. Since its discovery, the asymmetry in the opposite sides of the bipolar jet 
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has intrigued astronomers for many years. Previous analysis of the same STIS dataset suggested 
a nearly equa l mass-loss rate in th e opposite jets themselves but with asymmetry in the ambient 



environment (jMelnikov et al,l 12009 ). Through the demonstration of a sample fit with the X-wind 
models, we are able to interpret the asymmetry in velocity and intensity as a consequence of 
different mass loading on the opposite sides of the disk while the system adjusting itself to conserve 
linear momentum. Moreover, our analysis suggests the jet and wide-angle wind system as a whole 
is one entity, since observational interpretation through the optical lines alone would not recover 
the entirety of mass loss inside the hemisphere. 

In this work we demonstrate that the framework of X-winds can accommodate many gen- 
eral features of asymmetric disk-jet systems such as RW Aur A. Analysis such as one presented 
here may be applied to other similar systems. Further constraints of the disk-jet system through 
high-dispersion spectral line diagnosis associated with the jet and accretion activities are strongly 
encouraged. 

Facility: HST (STIS) 
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Table 1. Inferred Stellar Parameters 



Spec Type 
(1) 


logTeff 

(2) 


Band 
(3) 


m\ 
(4) 


(5) 


A v 
(6) 


(7) 


R* 

Rq 

(8) 


M* 
M 

(9) 


Reference 
(10) 


K4 


3.662 


J 


8.54 


0.18 


2.2 


0.36 


2.4 


0.85 


1 


K4 


3.662 








0.8 




2.0 


1.0 


2 


K3 


3.672 


J 


8.66 


0.49 


0.53 


0.05 


1.6 


1.3 


3 


K1/K3 


3.706/3.675 


Ic 


9.43 


~ 


0.39 


0.23 


1.7/2.2 


1.34/1.4 


4 



Note. — Description of columns 



Spectral type, with a typical uncertainty of one subclass; 



Effec tive temperature (in K) of the star for the spectral type in column (1) (jHartigan et al 



19941 1: 

The band used to calculate the stellar luminosity; the central wavelengths for J and Ic bands 
are 1.25 and 0.79 /im, respectively; 

The magnitude (in mag) of the star measured with the band in column (3); 
Veiling (defined by the ratio between total continuum flux and stellar flux) at the central 
wavelength of the band in column (3), derived from the veiling at 5700 A (ry) by assuming 
tj = 0.1 ry and that rj c is negligible; 

Visual extinction (in mag) used to calculate the extinction at the central wavelength of the 
band in column (3); 

Stellar luminosity (in solar luminosity) calculated from the magnitude in col umn (4), veiling 



i n col umn (5), extinction in column (6), and the bolometric corrections from lHartigan et al 



(|1994h : 

Stellar radius (in solar radius) derived from the effective temperature in column (2) and the 
luminosity in column (7); 

Stellar mass (in solar mass) inferred from the theoretical evolutionary tracks on the H-R 
diagram. 



Ref erenc es. - (1) lHartigan et al.l (|1995l ): (2) 
fcoOll ): (4) Iwhite fc Ghej fcoOll ~ 



Stout-Batalha et al 



(1200(1 ); (3) 



Woitas et al 
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Table 2. Radial Velocities in Literature and Resulting Velocity Ratios 



V Ta ,d 


min v rad 


max v vaxi 


Epochs 


^bluc 


Vred 


^blue/Ved 


Reference 


(km s" 1 ) 


(km s _1 ) 


(kms -1 ) 


(year) 


(kms" 1 ) 


(km s" 1 ) 






+14.0 ±4.9 


+7.3 


+22.4 


1981.9 - 1985.1 


-170.5 


+ 109.5 


1.56 


1 


+11.0 ±10.0 


-2.0 


+22.0 


1996.9 


-167.5 


+112.5 


1.49 


2 


+16.0 + 5.5 


+6.1 


+25.1 


1995.9 - 1999.9 


-172.5 


+107.5 


1.60 


3 


+16.0 + 1.9 






2003.2 


-172.5 


+107.5 


1.60 


4 


+23.5 






2000.9 


-180.0 


+100.0 


1.80 


5 



References. 



(120011 ): (4) IWhite fc Hillenbrand! (120041 '): (5) 



1) lHartmann et al. (11986T): (2) IStout-Batalha et all (120001 ): (3) IPetrov et al 



Woitas et al 



(2005). 



Table 3. Sample Fitting Set I 



i 


v w (blue)/v x 


v w (ved)/v x 




R x 


(°) 






(kms" 1 ) 


(AU) 


45 


3.208 -»• 3.609 


2.005 -»• 2.306 


70.534 


0.1779 


50 


3.208 -> 3.609 


2.005 -»• 2.306 


77.591 


0.1470 


55 


3.208 -)■ 3.609 


2.005 -»• 2.306 


86.953 


0.1170 


60 


3.208 -)■ 3.609 


2.005 -»• 2.306 


99.751 


0.0889 




Table 4. 


Sample Fitting 


Set II 




i 


v w (blue)/v x 


f w (red)/f x 




R x 


(°) 






(kms" 1 ) 


(AU) 



45 2.005 

50 2.005 

55 2.005 

60 2.005 



2.306 1.253 

2.306 1.253 

2.306 1.253 

2.306 1.253 



1.473 
1.473 
1.473 
1.473 



112.85 
124.15 
139.13 
159.60 



0.0695 
0.0574 
0.0457 
0.0347 
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Table 5. Adopted Parameters for the RW Aur General Set: (M*, #*) = (1M , 2R Q ) 



Properties 


Symbol 


Blue Jet 


Red Jet 


Unit 


Average terminal wind velocity 




310 


200 


kms -1 


Inclination angle 


i 


55 


55 


o 


Disk truncation radius 


R x 


0.117 


0.117 


AU 


Mass-loss rate 


M w 


3 x 10" 8 


5 x 10~ 8 


Moyr^ 1 


X-ray luminosity 


Lx 


6.0 x 10 31 


9.9 x 10 31 


ergs -1 


Mechanical heating coefficient 


a h 


1.0 x 10~ 3 


1.5 x 10~ 3 
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Fig. 1. — Narrow band images of the RW Aur A jet extracted from the HST/STIS archive. The jet 
axis is plotted as the abscissa, with the distance increasing along the position angle of 130°. The 
ordinate is the distance across the jet and is marked by positions of the seven slits separated by 
(X'07. The images are logarithmically scaled between 5 x 10~ 1 6 and 5 x 10~ 13 ergs -1 cm" 2 arcsec -2 . 



On each panel, the knots are identified and labeled as in lLopez-Martm et al.l ()2003l ). and the 
"stationary knot" on each side of the jet is labelled in parentheses by this work. 
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[OI]A6300 (Red), [0I]A6363 (Blue) 
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Fig. 2. — Position- Velocity diagrams of the central Of.'l in width of the RW Aur A jet derived from 
the HST/STIS spectra. Vr and Vb stand for the absolute values of the projected velocities of the 
red and the blue jet, respectively. Contours plotted start from 3a and increase by y/2a. The la of 
the lines are derived from the line-free regions next to the specific lines: 4.91 x 1CP 16 , 6.15 x 10~ 16 , 
5.15 x 10~ 16 , and 4.96 x 1(T 16 ergs" 1 cm" 2 arcsec" 2 for [O i] A6300, [O i] A6363, [N II] A6583, and 
[S n] A6731, respectively. 
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-1.0 -0.5 0.0 0.5 -1.0 -0.5 0.0 0.5 -1.0 -0.5 0.0 0.5 -1.0 -0.5 0.0 0.5 

log([SII]X6731/[OI]X6300) log([SII]A6731/[OI]X6300) log([SII]X6731/[OI]X6300) log([SII]X6731 /[OI]X6300) 



Fig. 3. — Optical line ratio diagrams for the RW Aur A jet derived from the HST/STIS spectra. 
The blue jet is defined here by the integrated intensities over velocity interval [—250, —20] kms -1 
and the red jet is defined by those within [+20, +200] kms -1 . On each panel, points labelled by the 
same color are sampled every / /05 unless the fractional error of either line ratio is above 50%. Left: 
Line ratios integrated over the central five slits at different distances along the jet. Right: Line 
ratios obtained at different positions from the central five slits without integration. The innermost 
/ /95 of the blue jet and / /65 of the red jet are shown for clarity. 
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20 40 20 40 



G3 (AU) 03 (AU) 

Fig. 4. — Profiles of the velocity (black contours, in unit of kms -1 ) and electron density (color 
map, in logarithmic unit of cm -3 ) of the X-wind solutions for the sample parameters listed in Table 
[5j Models adopted for the blue and red jets are shown on the left and right, respectively. 
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03 (AU) 03 (AU) 

Fig. 5. — Contours of electron fraction (black lines) are shown on top of the temperature profiles 
(color map, in unit of K), for the models listed in Table The left panel is for the blue jet and 
the right panel is for the red jet. 
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Fig. 6. — Images and spectra as observed and theoretical data are compared on top of each other 
in [O i] A6300 (for the red jet) and [O i] A6363 (for the blue jet). Left: Variations of the peak 
intensity and the full-width-half-maxium (FWHM) are plotted along the axes of jets, with their 
gaussian fits from the synthetic images (thick lines). The red jet is shown in red and the blue jet is 
shown in blue. Observed data points are shown in red filled circles for the red jet and in blue filled 
squares for the blue jet. Right: The position- velocity diagrams extracted from the observed (black 
contours) and synthetic (blue and red contours) spectra for the red and the blue jets in the upper 
panel. The Gaussian-fitted velocity centroids from the observed (blue and red filled symbols) and 
synthetic (blue and red thick lines) spectra are plotted in the lower panel, showing the distinct 
velocity ratios of about 1.6. 
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Fig. 7. — Images and spectra as observed and theoretical data are compared on top of each other 
in [N n] A6583. The meaning of the symbols is the same as that in Figure [6l 
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Fig. 8. — Images and spectra as observed and theoretical data are compared on top of each other 
in [S II] A6731. The meaning of the symbols is the same as that in Figure 
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Fig. 9. — Line ratios calculated from the blue jet and the red jet models (shaded areas) overlaid 
on the observed data (discrete symbols). The discrete symbols and shaded areas have the same 
spatial coverage in x ((X'12, roughly corresponding to the central three slits) and y (0'/95 for the 
blue jet and (X'65 for the red jet). 



